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Abstract: We study the capability of the MAGIC telescope to observe under moderate moonlight. TeV 7- 
ray signals from the Crab nebula were detected with the MAGIC telescope during periods when the Moon 
was above the horizon and during twilight. This was accomplished by increasing the trigger discriminator 
thresholds. No change is necessary in the high voltage settings since the camera PMTs were especially 
designed to avoid high currents. We characterize the telescope performance by studying the effect of the 
moonlight on the 7-ray detection efficiency and sensitivity, as well as on the energy threshold. 



Introduction 

Ground-based searches for very high energy 
(VHE) 7-ray emission from celestial objects are 
normally carried out by so-called imaging air 
Cherenkov telescopes (lACT) during clear, moon- 
less nights. If such a strict requirement is relaxed 
to allow observations under moderate moonlight or 
twilight, an increase of the duty cycle to 18% (from 
^ 1 000 to ^ 1 500 hours of observations per year) is 
possible. VHE observations under moonlight have 
been tested in the past [8, 4, 6] but with solutions 
that were too expensive, time consuming and not 
efficient in terms of energy threshold and sensitiv- 
ity. The MAGIC lACT [7] has been designed to 
carry out observations also during moderate moon- 
light. This places MAGIC in a prominent posi- 
tion, in particular for the study of variable sources 
as well as in multi-wavelength campaigns together 
with other instruments. In this paper we describe 
the technical innovations and analysis changes that 
allow observations in the presence of the Moon 
with MAGIC. 

Traditionally, PMTs are operated with amplifica- 
tion gains around 10^ — 10^ which, under moon- 
light, generate continuous (direct) currents (DCs) 
that can damage the last dynode, resulting in rapid 
ageing of the PMT. MAGIC PTMs run at a gain of 



about 3 X 10^ [2] . In order to also detect single pho- 
toelectrons (phe) the PMT signal is fed to an AC- 
coupled fast, low noise preamplifier to raise the 
combined gain to about 10^. The PMT analog sig- 
nal is transmitted over an optical fiber, converted 
into an electrical pulse and split into two branches, 
one of which enters a discriminator (DT), which 
issues a digital signal whenever the pulse exceeds 
a given threshold. 

The increase of the background light due to the 
presence of the Moon depends on various factors 
including Moon and source zenith angles. Moon 
phase, angular distance to the Moon and atmo- 
spheric conditions. At around 25° away from the 
Moon the direct scattered moonlight approaches 
a constant level below that of the night sky light 
background. Figure 1 shows the dependence of the 
trigger rate (for a four neighboring pixels config- 
uration) on the DT settings, for different illumi- 
nation conditions (which produce different anode 
currents in the camera). The light of night sky 
(LONS) is responsible for the steep increase at low 
DT values (at ^30 a.u. in the case of dark obser- 
vations). At higher DT values the rate is caused 
mostly by Cherenkov showers. The telescope op- 
erates at the minimum possible DT for which the 
contribution of accidental triggers is negligible. 
For extragalactic regions (DC=1//A) the DTs are 
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Figure 1 : Trigger rate as a function of the discrim- 
inator threshold for four neighboring pixels con- 
figuration and different camera illuminations. The 
shaded area shows the range used for MAGIC reg- 
ular observations (dark and under moonlight). The 
dashed line shows the linear regime. The upper 
axis shows the corresponding energy threshold (af- 
ter image cleaning) for observations at zenith an- 
gles between 20° and 30° as deduced in Section 3. 



generally set to 30 a.u., which corresponds to a 
pulse charge of 8-10 phe. Higher DT values are 
needed to keep the trigger rate below the limit of 
the DAQ system (500 Hz) for observations during 
twilight and moonlight. We restrict MAGIC ob- 
servations to a maximum DC of 8 /j,A. This per- 
mits observations in the presence of the Moon un- 
til (since) 3-4 days before (after) full Moon, for an 
angular distance to the Moon greater than 50° [2]. 



- Dark, DC 1 .0 fiA, DT 35 a.u. 
Weak moonlight, DC 2.4 fiA, DT47 a.u. 
Strang moonlight. DC 4.5 tiA, DT 60 a.u. 
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Figure 2: Distributions of SIZE before analysis 
cuts for three Crab nebula samples acquired un- 
der different light conditions. The histograms have 
been normalized to a common observation time. 
Note that the distributions are completely domi- 
nated by hadronic events (~ 99%). 



based on the trigger and after-cleaning rates were 
applied in order to remove bad runs. The shower 
images were parameterized using the Hillas param- 
eters [5] SIZE, WIDTH, LENGTH, DIST, CONC 
and ALPHA, combined (except for the latter) into 
an adimensional variable (HADRONNESS) for 
7/hadron separation by means of a Random For- 
est classification algorithm [3]. The signal region 
was defined by the cuts HADRONNESS<0.15 and 
ALPHA<8°. 



Results 



Observations and Data Analysis 

To characterize the response of the telescope un- 
der moonlight, we observed the Crab nebula at 
different light conditions (including dark obser- 
vations used as reference) between January and 
March 2006, in the ON/OFF mode. Two data 
sets, one with zenith angle between 20° and 30°, 
and a second one between 30° and 40°, were ac- 
quired and analyzed separately. Depending on 
the different moonlight levels, the resulting an- 
ode currents ranged between 1 and 6 /iA. Corre- 
spondingly, the DT was varied between 35 and 
65 units. The acquired data were processed by 
the standard MAGIC analysis chain [1]. The im- 
ages were cleaned using absolute tail and bound- 
ary cuts at 10 and 5 phe, respectively. Quality cuts 



As we increase the DT levels to counteract acci- 
dental triggers, one depletes the SIZE distribution 
of shower candidates, as expected, mostly at low 
values. However, we find that a substantial num- 
ber of showers with SIZE up to 10'' phe, i.e. those 
well above the trigger level (which is around 50 
phe), are also suppressed (see Figure 2). On the 
other hand, the LENGTH, WIDTH and CONC dis- 
tributions above 200 phe do not show significant 
differences (See Figure 3) and therefore we do not 
expect 7/hadron separation to degrade due to the 
presence of the Moon. Below 200 phe the distribu- 
tions are distorted by the different trigger thresh- 
old. 

The supressed high-SIZE showers are those at high 
DIST (see Figure 3c), corresponding to a high im- 
pact parameter, where a significant fraction of the 
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Figure 4: Relative 7-ray detection efficiency 
(green triangles, left axis) and sensitivity (red cir- 
cles, right axis) as a function of DT (SIZE> 
400 phe), for zenith angle bins [20°, 30°] (filled 
markers) and [30°, 40°] (empty markers) measured 
from Crab nebula observations. 



Figure 3: Distributions of LENGTH (a), WIDTH 
(b) and DIST for all recorded events (c) and for 
images fully contained in the inner camera (d) for 
SIZE>400 phe. Three Crab nebula samples ac- 
quired under different moonlight conditions. The 
histograms are normalized to a unit area in (a) and 
(b) and to a common observation time in (c) and 
(d). 



light falls outside the trigger area ( 1 ° radius around 
the camera center). In some cases, the fraction 
of the shower image contained in the trigger area 
does not exceed the increased threshold for at least 
4 neighboring pixels, as required for a trigger. A 
confirmation of this hypothesis was obtained by 
two different tests. First, the DIST distributions 
for events fully contained in the trigger region were 
compared (Figure 3d). In such a case we find simi- 
lar distributions. The second test was performed by 
observing Crab in dark conditions (DC^l.l /iA), 
but with increased DTs. In this case we found 
similar inefficiencies as those shown in Figure 3c). 
Therefore we can conclude that the change of the 
DIST distribution is not related to the mean DC 
current (i.e. with the camera illumination) but only 
to the DT level. 

These results show that moonlight does not distort 
the images from Cherenkov showers and therefore 
the analysis based on the Hillas parameters does 
not have to be adapted for data acquired under 
moonlight, and in particular the 7/hadron separa- 
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Figure 5: Energy threshold after image cleaning as 
a function of DT obtained from MC simulated 7- 
ray events (for zenith angle between 20° and 30°). 
The top axis shows the typical mean DC for a cho- 
sen DT value. 



tion power is not reduced for this kind of obser- 
vations. In addition, the differences that we find 
in the event rates and the DIST distributions are 
exclusively due to the fact that the DTs were in- 
creased to keep a low rate of accidental events, to- 
gether with the fact that the trigger area does not 
span the whole camera. With a faster DAQ sys- 
tem one could have a fix DT level and deal with 
the extra noisy events introduced by the moonlight 
during the offline analysis. 

For this study the dependence of the telescope re- 
sponse on the level of Moon illumination will be 
parameterized as a function of DT. As is shown 
in Figure 2, there is a reduction in the collection 
area over a wide range of energies. We use Crab 
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nebula observations to parametrize the efficiency 
of detecting 7-rays, for every DT and SIZE, rel- 
ative to the values for dark observations. These 
values are used during the off-line analysis to- 
gether with the MC simulation with standard DTs 
to calculate the correct collection areas. For this, 
the observations of the Crab nebula are divided 
into different samples according to the observa- 
tion date, SIZE and DT values. For each of the 
samples we get a measurement of the 7-ray rate 
(i?) and sensivity (.s). We find that, for a given 
SIZE range, the dependences of R/Rq and s/sq 
with the DT are well described by linear functions 
(i?/i?o, s/so = (1 - Sr,s [DT - DTo))), where 
Ro, So and DTq are the rate, sensitivity and DT 
values, respectively, for the dark observation case. 
Sn^a are let free during the fit and are referred to 
as efficiency/sensitivity loss rate, respectively. The 
results for SIZE>400 phe and the two considered 
zenith angle samples ([20°, 30°] and [30°, 40°]) are 
shown in Figure 4. The fit parameters obtained for 
both zenith angles are compatible within statistical 
errors. This allows us to perform a combined fit 
for both samples. For SIZE>400 phe we obtain 
the following expression for the 7-ray efficiency 
and sensitivity: 

R/Ro = 1 - (1.41 ± 0.32) X IQ-^pT - DTq) 
s/so = 1-F(6.3±1.6) X 10-^(£)T-£)To) 

In order to understand the dependence of the 7- 
ray detection efficiency on the energy we have 
performed the same study for four bins of SIZE, 
namely [200,400], [400,800], [800,1600] and 
[1600,6400] phe, which roughly correspond to the 
energy ranges 1 150,300], [300,600], [600-1000] 
and > 1000 GeV, respectively, for low zenith an- 
gle. Up to SIZE=3000 phe we find a hnear depen- 
dence that can be parameterized by: 

Sr = (2.24 ± 0.13) X 10-2 

-(7.2 ± 1.2) X 10"*^ SIZE [phe] 
Sa = (1.63 ± 0.14) X 10^2 

-(7.4 ± 1.8) X 10"*^ SIZE [phe] 

We tested the dependences of these results on dif- 
ferent HADRONNESS and ALPHA cuts, obtain- 
ing in all cases similar results. 

Finally, it is important to understand the influence 
of the moonlight on the energy threshold. We de- 



fine the energy threshold as the peak of the energy 
distribution of all events after image cleaning and 
before analysis cuts. The dependence of the energy 
threshold is well described by the foUowing hnear 
function (see Figure 5): 

E 

— ^ = (69.3±0.4)-|- (1.06 ±0.03) (DT-DTq) 
GeV 

This increase is relatively marginal, and it has to 
be noted again that it is due to the increase of the 
DTs, and hence only indirectly to the increase in 
the camera illumination. 

We would like to thank the lAC for the excellent 
working conditions at the Observatorio del Roque 
de los Muchachos in La Palma. 
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